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Abstract. We study some global aspects of the bifurcation of an equivariant family of volume-contracting 
vector fields on the three-dimensional sphere. When part of the symmetry is broken, the vector fields 
exhibit Bykov cycles. Close to the symmetry, we investigate the mechanism of the emergence of het- 
eroclinic tangencies coexisting with transverse connections. We find persistent suspended horseshoes 
accompanied by attracting periodic trajectories with long periods. 



1. Introduction 

Heteroclinic cycles and networks associated to equilibria, periodic solutions and chaotic sets may be 
responsible for intermittent dynamics in nonlinear systems. Heteroclinic cycles may also be seen as 
the skeleton for the understanding of complicated switching between physical states - see Field [TU], 
Golubitsky and Stewart [TT] and Melbourne et al [TO] . 

The homoclinic cycle associated to a saddle-focus [T3] provides one of the main examples for the 
occurrence of chaos involving suspended hyperbolic horseshoes and strange attractors; the complexity 
of the dynamics near these cycles has been discovered by the pionner L. P. Shilnikov (271 [28l [SO] ■ The 
simplest heteroclinic cycles between two saddle-foci of different Morse indices where one heteroclinic 
connection is structurally stable and the other is not have been first studied by Bykov fSJ and are thus 
called Bykov cycles. Recently there has been a renewal of interest of this type of heteroclinic bifurcation 
in different contexts - see [TU [T71 [M] and references therein. We also refer Lamb et al [15] who have 
studied Bykov cycles in the context of reversible systems. 

Explicit examples of vector fields for which such cycles may be found are reported in Aguiar et al [3] 
and Rodrigues and Labouriau [25 . These examples start with a differential equation with symmetry, 
X = fo{x) whose fiow has a globally attracting three-dimensional sphere, containing an asymptotically 
stable heteroclinic network with two saddle-foci. When part of the symmetry is destroyed by a small non- 
equivariant perturbation, it may be shown by the Melnikov method that the two-dimensional invariant 
manifolds intersect transversely. When some symmetry remains, the connection of the one-dimensonal 
manifolds is preserved, giving rise to Bykov cycles forming a network. 

The main goal of this article is to describe and characterize the transition from the dynamics of the 
flow of the fully symmetric system x = fo{x) and the perturbed system x = f\{x), for small A 7^ 0. For 
X = fo{x), there is a heteroclinic network whose basin of attraction has positive Lebesgue measure. 
When A 7^ 0, the intersection of the invariant manifolds is transverse giving rise to a network E* that 
cannot be removed by any small smooth perturbation. The transverse intersection implies that the set 
of all trajectories that lie for all time in a small neighbourhood of Y,* has a locally-maximal hyperbolic 
set, admiting a complete description in terms of symbolic dynamics [53]. Labouriau and Rodrigues [T7] 
proved that for the perturbed system x = f\{x), the fiow contains a Bykov network E* and uniformly 
hyperbolic horseshoes accumulating on it. 

Suppose the fully symmetric network E*' is asymptotically stable and let V'~' be a neighbourhood of 
E'' whose closure is compact and positively fiow- invariant; hence it contains the w-limit sets of all its 
trajectories. The union of these limit sets is a maximal invariant set in V^. For x — foix), this union 
is simply the network Eq. For symmetry-breaking perturbations of /o it contains, but does not coincide 
with, a nonwandering set r2A(E*) of trajectories that remain close to E*, the suspension of horseshoes 
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accumulating on T,*. The goal of this article is to investigate the larger limit set that contains nontrivial 
hyperbolic subsets and attracting limit cycles with long periods in V'^. This is what Gonchenko et al 
[12) call a strange attractor: an attracting limit set containing nontrivial hyperbolic subsets as well as 
attracting periodic solutions of extremely long periods. 

When A 0, the horseshoes in £7^(2*) lose hyperbolicity giving rise to heteroclinic tangencies with 
infinitely many sinks nearby. A classical problem in this context is the study of heteroclinic bifurcations 
that lead to the birth of stable periodic sinks - see Afraimovich and Shilnikov [T] and Newhouse [HJ [22] • 
When we deal with a heteroclinic tangency of the invariant manifolds, the description of all solutions 
that lie near the cycle for all time becomes more difficult. The problem of a complete description is 
unsolvable: the source of the difficulty is that arbitrarily small perturbations of any differential equation 
with a quadratic homo/heteroclinic tangency (the simplest situation) may lead to the creation of new 
tangencies of higher order, and to the birth of degenerate periodic orbits — Gonchenko [T5] . 

Large-scale invariant sets of planar Poincare maps vary discontinuously in size under small pertur- 
bations. Global bifurcations of observable sets, such as the emergence of attractors or metamorphoses 
of their basin boundaries, are easily detected numerically and regularly described. However, in the ex- 
ample described in [25], the global bifurcation from a neighbourhood of to V'^ is still a big mistery. 
The present paper contributes to a better understanding of the transition between uniform hyperbolicity 
(Smale horseshoes with infinitely many slabs) and the emergence of heteroclinic tangencies in a dissipative 
system close the symmetry. 

Pramework of the paper. This paper is organised as follows. In section [3] we state our main result 
and review some of our recent results related to the object of study, after some basic definitions given in 
section [2] The coordinates and other notation used in the rest of the article are presented in section [H 
where we also obtain a geometrical description of the way the flow transforms a curve of initial conditions 
lying across the stable manifold of an equilibrium. In section [5j we investigate the limit set that contains 
nontrivial hyperbolic subsets and we explain how the horseshoes in ^Ia(S*) lose hyperbolicity, as A — ^ 
0. The first obstacle towards hyperbolicity is the emergence of tangencies and the existence of thick 
suspended Cantor sets near the network. In section [HI we prove that there is a sequence of parameter 
values Xi accumulating on such that the flow of f\- has heteroclinic tangencies and thus infinitely many 
attracting periodic trajectories. We include in section [7] a short conclusion about the results. 

2. Preliminaries 

Let / be a C"'^ vector field on R" with flow given by the unique solution x(t) — (p{t,xo) € R" of 

1 = f{x) and a;(0) = xq. Given two equilibria p and an m-dimensional heteroclinic connection from p 
to denoted [p — > g], is an m-dimcnsional connected flow-invariant manifold contained in {p)r\W'^ {q) . 
There may be more than one trajectory connecting p and q. 

Let S —{pj : j G {!,..., k}} be a flnite ordered set of mutually disjoint invariant saddles. Following 
Field |10[ . we say that there is a heteroclinic cycle associated to S if 

Vj e {1, . . . , k}, W"{p/) n W'{pj+i) ^ (mod k). 

Sometimes we refer to the equilibria defining the heteroclinic cycle as nodes. A heteroclinic network is 
a finite connected union of heteroclinic cycles. Throughout this article, all nodes will be hyperbolic; the 
dimension of the local unstable manifold of an equilibria p will be called the Morse index of p. 

In a three-dimensional manifold, a Bykov cycle is a heteroclinic cycle associated to two hyperbolic 
saddle-foci with different Morse indices, in which the one-dimensional manifolds coincide and the two- 
dimensional invariant manifolds have a transverse intersection. It arises as a bifurcation of codimension 

2 and it is also called by T -point. 

Let II denote a measure on a smooth manifold AI locally equivalent to the Lebesgue measure on charts. 
Given a; S M, let 

ij{x) = fl W{t,x):t>T} 

T>0 

denote the oj-limit set of the solution through x. If X C M is a compact and flow-invariant subset, 
we let B{X) — {x £ M : uj{x) C X} denote the basin of attraction of X. A compact invariant subset 
X of M is a Milnor attractor if ii{B{X)) > and for any proper compact invariant subset Y of X, 
H{B{X)\B(Y)) > 0. 
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Figure 1 . There are two different possibilities for the geometry of the flow around a Bykov 
cycle depending on the direction trajectories turn around the heteroclinic connection [v ^ w]. 
We assume here that all trajectories turn in the same direction near v and near w as drawn on 
the left (a) . When the endpoints of the trajectory are joined, the closed curve is linked to the 
cycle. The case where trajectories turn in opposite directions around the connection is shown 
in (b). In this case, when the endpoints of the trajectory are joined, the closed curve may not 
be linked to the cycle. 



3.1. The organising centre. The starting point of the analysis is a differential equation i = foix) on 
the unit sphere — {X = {xi,X2,X'i,Xi) e R"* : \\X\\ = 1} where /q-.S^^ TS^ is a vector field 
with the following properties: 

(PI) The vector field /o is equivariant under the action of Z2 © Z2 on S'^ induced by the two linear 
maps on R*: 

7i(xi,X2,a;3,a;4) = (-xi, -X2, ^3, 2:4) and 72(2^1, 2:2, a;3, ^4) = {xi, X2, -X3, x^). 

(P2) The set Fix{7i2 © Z2) = {x S S'^ : 71a; = 72a; = x} consists of two equihbria v = (0, 0, 0, 1) and 
w = (0, 0, 0, —1) that are hyperbolic saddle-foci, where: 

• the eigenvalues of d/o(v) are — Cv ± a^i and with av 7^ 0, Cv > -Ev > 

• the eigenvalues of dfo{w) are E-^ ± aw* and — Cw with aw 7^ 0, Cw > -Ew > 0. 

(P3) The fiow- invariant circle Fix{{'yi)) = {a; G S'^ : 71a; = a;} consists of the two equilibria v and w, 
a source and a sink, respectively, and two heteroclinic trajectories from v to w that we denote 
by [v — > w] . 

(P4) The /o-invariant sphere Fix{{-f2)) = {a; G : 72a; — x} consists of the two equilibria v and w, 
and a two-dimensional heteroclinic connection from w to v. Together with the connections in 
(F|2|) this forms a heteroclinic network that we denote by S". 

(P5) For sufficiently small open neighbourhoods V and of v and w, respectively, given any trajectory 
(f going once from V to W, if one joins the starting point of (p in dV to the end point in dW by 
a line segment, one obtains a closed curve that is linked to (figure [T|). 
Condition (F|3]) means that the curve ip and the cycle E° cannot be separated by an isotopy. This 
property is persistent under perturbations: if it holds for the organising centre /o, then it is still valid 
for vector fields near it, as long as the heteroclinic connection remains. An explicit example of a family 
of differential equations where this assumption is valid is constructed in Rodrigues and Labouriau |25| . 

3.2. The heteroclinic network of the organising centre. The heteroclinic connections in the net- 
work are contained in fixed point subspaces satisfying the hypothesis (HI) of Krupa and Melbourne 
|16) . Since the inequality CvCw > E^E^ holds, the Krupa and Melbourne stability criterion [TB] may 
be applied to E° and we have: 

Proposition 1. Under conditions (in})-(I^ the heteroclinic network E*' is asymptotically stable. 
In particular, we obtain: 
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Corollary 2. The basin of attraction of the heteroclinic network S*^ whose vector field satisfies (I^^-flJ^ 
has positive Lebesgue measure. 

Proposition [1] implies that there exists an open neighbourhood V'^ of the network I]° such that every 
trajectory starting in V'^ remains in it for all positive time and is forward asymptotic to the network. 
The neighbourhood may be taken to have its boundary transverse to the vector field /q. 

The fixed point hyperplane defined by Fix{{'j2)) — {(a;i, X2, 2:3, X4) e S'^ : 0:3 = 0} divides S'^ in 
two flow- invariant connected components, preventing jumps between the two cycles in TP. Due to the 
(72)-symmetry, trajectories whose initial condition lies outside the invariant subspaces will approach one 
of the cycles in positive time. Successive visits to both cycles require breaking this symmetry. 

3.3. Bykov cycles. Since v and w are hyperbolic equilibria, then any vector field close to it in the 
topology still has two equilibria v and w with eigenvalues satisfying (F|2]). The dimensions of the 
local stable and unstable manifolds of v and w do not change. If we retain the symmetry 71, then the 
one-dimensional connection of (FO remains, as it takes place in the fiow-invariant circle Fix{{'yi)), but 
generically when the symmetry 72 is broken, the two dimensional heteroclinic connection is destroyed, 
since the fixed point subset _Fia;(Z2((72))) is no longer flow- invariant. Generically, the invariant two- 
dimensional manifolds meet transversely at two trajectories, for small (7i}-equivariant perturbations of 
the vector field. We denote these heteroclinic connections by [w v] This gives rise to a network E* 
consisting of four copies of the simplest heteroclinic cycle between two saddle-foci of different Morse 
indices, where one heteroclinic connection is structurally stable and the other is not. This cycle is called 
a Bykov cycle. Note that the networks and E* are not of the same nature. 

Transversality ensures that the neighbourhood is still positively invariant for vector fields close 
to fo and contains the network E*. Since the closure of V'^ is compact and positively invariant it contains 
the w-limit sets of all its trajectories. The union of these limit sets is a maximal invariant set in V'^ . For 
/o, this is the cycle E'^, by Proposition [TJ whereas for symmetry-breaking perturbations of fo it contains 
E* but does not coincide with it. From now on, our aim is to obtain information on this invariant set. 

A systematic study of the dynamics in a neighbourhood of the Bykov cycles in E* was carried out 
in Aguiar et aJ |2], Labouriau and Rodrigues [17] and Rodrigues [24]; we proceed to review these local 
results. In the next section we will discuss some global aspects of the dynamics. 

Given a Bykov cycle F involving v and w, let Uv and ?7w C be disjoint neighbourhoods of these 
points. Let Sp and Sq be local cross-sections of f\ at two points p and q in the connections [v — w] and 
[w — ^ v], respectively, with p,q ^ Uv U U^. Saturating the cross-sections by the fiow, one obtains two 
flow-invariant tubes joining and C/w containing the connections in their interior. We call the union of 
these tubes, and [/^ a tubular neighbourhood of the Bykov cycle. 

With these conventions we have: 

Theorem 3 (Labouriau and Rodrigues [T7|, 2012). // a vector field fo satisfies (in\)~(lWj), then the 
following properties are satisfied generically by vector fields in an open neighbourhood of fo in the space 
of {'yi)-equivariant vector fields of class on S'^ : 

(1) there is a heteroclinic network E* consisting of four Bykov cycles involving two equilibria v and 
w, two heteroclinic connections [v w] and two heteroclinic connections [w — >■ v]; 

(2) the only heteroclinic connections from to w are those in the Bykov cycles and there are no 
homoclinic connections; 

(3) any tubular neighbourhood of a Bykov cycle F in E* contains points not lying on F whose 
trajectories remain in for all time; 

(4) any tubular neighbourhood of a Bykov cycle F in E* contains at least one n-pulse heteroclinic 
connection [w v]; 

(5) given a cross-section Sq C at a point g m [w — S- v] , there exist sets of points such that the 
dynamics of the first return to Sq is uniformly hyperbolic and conjugate to a full shift over a finite 
number of symbols. These sets accumulate on the cycle. 

Notice that assertion ^ of Theorem [3] implies the existence of a bigger network: beyond the original 
transverse connections [w v], there exist infinitely many subsidiary heteroclinic connections turning 
around the original Bykov cycle. We will restrict our study to one cycle. 

It is a folklore result that a hyperbolic invariant set of a C^-diffeomorphism has zero Lebesgue measure 
- see Bowen However, since the authors of [T7] worked in the category, this chain of horseshoes 
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might have positive Lebesgue measure as the "fat Bowen horseshoe" described in [6 . Rodrigues pi] 
proved that this is not the case: 

Theorem 4 (Rodrigues |24) . 2013). Let be a tubular neighbourhood of one of the Bykov cycles T of 
Theorem\^ Then in any cross-section Sq C at a point g m [w — > v] the set of initial conditions in 
Sq n that do not leave for all time has zero Lebesgue measure. 

It follows from Theorem[4]that the shift dynamics does not trap most trajectories in the neighbourhood 
of the cycle. In particular, the cycle cannot be Lyapunov stable (and therefore cannot be asymptotically 
stable). 

One astonishing property of the heteroclinic network S* is the possibility of shadowing S* by the 
property called switching: any infinite sequence of pseudo-orbits defined by admissible heteroclinic 
connections can be shadowed, as we proceed to define. 

A path on S* is an infinite sequence s^ — {cj)jeT<! of heteroclinic connections Cj = [Aj Bj] in S* 
such that Aj,Bj € {v,w} and Bj = Aj+i. Let Uv,Uvr C V'^ be neighbourhoods of v and w. For each 
heteroclinic connection in E*, consider a point p on it and a small neighbourhood Up C V''^ of p. We 
assume that all these neighbourhoods are pairwise disjoint. 

The trajectory (p{t, q) oi i = f\{x), follows the path s*"' = (cj)jgN within a given set of neighbourhoods 
as above, if there exist two monotonically increasing sequences of times (fi)igN and (zi)igN such that for 
all i G N, we have ti < Zi < ti+i and: 

• Lp{t, q) C V° for all t > 0; 

• ip{ti,q) e Ua, and (p{zi,q) eUp, p e Q and 

• for all t € {zi, Zi^i), ip{t, q) does not visit the neighbourhood of any other node except 

There is switching near E* if for each path there is a trajectory that follows it within every set of 
neighbourhoods as above. 

Theorem 5 (Aguiar et aJ 2005). There is switching on the network S*. 

The solutions that realise switching lie for all positive time in the union of tubular neighbourhoods 
of all cycles F C S*. Hence we may adapt the proof of Theorem U] to obtain: 

Corollary 6. The switching of Theorem\^ is realised by a set of initial conditions with zero Lebesgue 
measure. 

3.4. Non-hyperbolic dynamics. When the symmetry 72 is broken, the dynamics changes dramatically, 
as can be seen in the next result: 

Theorem 7. // a vector field fo satisfies then in any open neighbourhood of fo in the space 

of {ji)-equivariant vector fields of class on S^, there are vector fields /* whose flow has a heteroclinic 
tangency between VF"(w) and ^{v) in V'^ . 

Theorem [7] is proved in section [51 The vector fields will be obtained in a generic one-parameter 
unfolding f\ of fo, for which we will find a sequence of Xi converging to zero such that the flow of f\- 
has the required property. When Xi — >■ 0, these tangencies accumulate on the transverse connections. 
Persistent tangencies in a dissipative diffeomorphism are related to the coexistence of infinitely many sinks 
and sources |21U22| . Moreover, any parametrised family of diffeomorphisms going through a heteroclinic 
tangency associated to a dissipative cycle must contain a sequence of Henon-like families ^ . Hence, for 
A w 0, return maps to appropriate domains close to the tangency are conjugate to Henon-like maps and 
thus: 

Corollary 8. Suppose a vector field fo satisfies Then in the space of {ji) -equivariant vector 

fields of class on S'^, there is a set C of vector fields accumulating on fo such that all f\ £ C possess 
infinitely many strange (coexisting) attractors in V'^ . 

Proof. In Ovsyannikov and Shilnikov |23j . it is shown that there are small perturbations of fo with a 
periodic solution as close as desired to the cycle whose stable and unstable manifolds are tangent. By 
Colli [9], the result follows. Note that although the statement in [9] asks for C°° perturbations, in [5] the 
coexistence of strange attractors only requires . □ 
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ln(w) Out(v,+) 



Figure 2. Cylindrical neighbourhoods of the saddle-foci w (a) and v (b). 

4. Local dynamics near the network 

In order to obtain results on (71 )-equi variant vector fields close to /o we study the bifurcation 
of a generic one-parameter family of differential equations x ~ f\{x) on the unit sphere = {X = 
(xi,a:2,a;3,X4) G R'' : \\X\\ = 1}. The unfolding fx:S^^ TS^ of fo is a family of vector fields with 
the following properties: 

(P6) For each A the vector field f\ is (7i)-equivariant. 

(P7) There are two equilibria v and w satisfying (F(5]) and (F^. 

(P8) For A ^ 0, the local two-dimensional manifolds VF"(w) and W^iv) intersect transversely at two 
trajectories. Together with the connections in (I^ this forms a Bykov heteroclinic network that 
we denote by T,^. 

In order to describe the dynamics around the Bykov cycles, we start by introducing local coordinates 
near the saddle-foci v and w and we define some terminology that will be used in the rest of the paper. 
Since by assumption (F(2]) we have Cv 7^ and Cw 7^ E^, then by Samovol's Theorem [53], the vector 
field fx is C^~conjugate to its linear part around each saddle-focus — see also Homburg and Sandstede 
[T5] (section 3.1). In cylindrical coordinates {p,9,z) the linearization at v takes the form 

p = —CvP = 1 z = EvZ 

and around w it is given by: 

p = E^p 9 = 1 z = —Cv,z. 

We consider cylindrical neighbourhoods V and W in S'^ of v and w, respectively, of radius e > 
and height 2e — see figured After a linear rescaling of the local variables, we may take e = 1. Their 
boundaries consist of three components: the cylinder wall parametrized by a; G R (mod 2tt) and |?/| < 1 
with the usual cover (x, y) ^ (1, x, y) = {p, 0, z) and two discs (top and bottom of the cylinder). We take 
polar coverings of these discs {r,(p) 1— {r,tp,±l) = {p,0,z) whith < r < 1 and G R (mod 2tt) and 
use the following terminology, as in figure [21 

• /n(v), the cylinder wall of V, consists of points that go inside V in positive time; 

• Out{v), the top and bottom of V, consists of points that go inside V in negative time; 

• /n(w), the top and bottom of W, consists of points that go inside W in positive time; 

• Out{w), the cylinder wall of W, consists of points that go inside W in negative time. 

The flow is transverse to these sets and moreover the boundaries of V and of W may be written as the 
closures of the disjoint unions /n(v) U Out{v) and /n(w) U Out{^v), respectively. The trajectories of all 
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Figure 3. For A close to zero, we are assuming that VK°(v) intersects tiie wall Out{w) of 
the cylinder W on a, closed curve - this is the expected unfolding from the coincidence of the 
invariant manifolds of the equilibria. 



points {x,y) in In{\-)\W'' (v), leave V at Out(v) at 



(1) = |^|y|^.,__Ji^+^j = (^,0) where ^ ^ > 1 

Similarly, points (r, (p) in /n(w)\M^*(w), leave W at Out{w) at 

In r 5„ ^ ' ^ ■ - 



(2) $w(r, (p) = (^-— + if, r'- j - (a:, y) where 6^ = — > I . 

We will denote by Wl^^{v) the portion of W"{v) that goes from v up to In{w) not intersecting the 
interior of W. Similarly, Wf^^{v) is the portion of W^{v) outside W that goes directly from Out{w) into 
V, and W^^{'w) and Wf^^{'w) connect w to dV, not intersecting the interior of V. 

The flow sends points in Out{v) near Wj"^(v) into In{w) along the connection [v — !■ w]. We will 
assume that this map is the identity, this is compatible with hypothesis (FO; nevertheless all the 

results follow if ^v-s-w is either a uniform contraction or a uniform expansion. We make the convention 
that one of the connections [v w] links points with i/ > in 1/ to points with y > in W. There is 
also a well defined transition map : Out{w) — > /n(v) that will be discussed later. 

By (F[S|), the manifolds M^"(w) and ly (v) intersect transversely for A 7^ 0. For A close to zero, we 
are assuming that Vl^^(v) intersects the wall Out{w) of the cylinder on a closed curve represented in 
figure [3] by an ellipse — this is the expected unfolding from the coincidence of the invariant manifolds of 
the equilibria. 
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From the geometrical behaviour of the local transition maps ^ and we need some definitions: 
a segment (5 on In(y) is a smooth regular parametrized curve /3 : [0,1) — In{'v) that meets Wf^^lv) 
transversely at the point /?(!) only and such that, writing /3(s) = {x{s), y{s)), both x and y are monotonic 
functions of s - see figure [H] (a). 

A spiral on a disc D around a point p G 13 is a curve a : [0, 1) — > -D satisfying lim a{s) = p and 

such that if a{s) — (r(s), 0{s)) is its expressions in polar coordinates around p then the maps r and 9 are 
monotonic, and lim^^^- = +oo. 

Consider a cylinder C parametrized by a covering G R x [a,b], with a < 6 G R where 6' is 

periodic. A helix on the cylinder C accumulating on the circle h = Hq is a curve 7 : [0, 1) — 5- C such that 
its coordinates {9(s), h(s)) satisfy lim^^j^- h{s) = Hq, lim^^i- = +00 and the maps 9 and h are 

monotonic. 

Using these definitions and the expressions ((T|) and ^ for <i>v and <i>w we get: 

Proposition 9 (Aguiar et al 4J, 2010). When f/f^j holds, then a segment on In{v) is mapped by $v 
into a spiral on Out{'v) around W;"g(v) n Oui(v). This spiral is mapped by \l/v^.w into another spiral 
around Wf^^{'w) D In{v), which is mapped by $w into a helix on Out(w) accumulating on the circle 

OMt(w) n W"(w). 

5. Hyperbolicity 

In this section we show that the hyperbolicity of Theorem [3] only holds when we restrict our attention 
to trajectories that remain near the cycles in the network. The construction also indicates how the 
geometrical content of Theorems [3l |4] and [5] is obtained. 

Let (P.^,0) and (P.^,0) with < < < 27r be the coordinates of the two points in W;"c(w) n 
W;*^(v)nOui(w) where the connections [w — > v] meet Ouf(w), as in figure[31 Analogously, let (P^, 0) and 
(P^, 0) be the coordinates of the two corresponding points in Wj"^(w) n Wf^^{v) fl In(v) where [w — > v] 
meets In{v), with the convention that (P^, 0) and (P^, 0) are on the same trajectory for j = 1, 2. 

For small A > we may write W'(v) n Out{w) as the graph of a smooth function y — g{x), with 
g(Pi,) = 0, j = 1,2. Similarly, Vl^"(w) n In{v) is the graph of a smooth function y = h{x), with 
g{Pl) ^ 0, j ~ 1, 2. For definiteness, we number the points in the connections to have g'{P^) > 0. Hence 
g'iPl) < 0, and h'{P^) < 0, /i'(p2) > 0. 

Proposition 10. For the first hit map rj = <i>w o ^'v^w ° $v : In{v) — > Out{w), we have: 

i) any horizontal line segment [a, b] x {yo} C In{v) is mapped by 77 into a horizontal line segment 
[c,d] X {2/0} C Out{w), with S > 1; 

ii) any vertical line segment {xq} x [0, yo] C /n(v) is mapped by rj into a helix accumulating on the 
circle Out(w) n Wil^{w); 

iii) given xq, there are positive constants a < G R and a sequence of intervals 

1-n = {a;o} X [e^^^'^Z-'^e", e~2n7r/Kg6j ^^^j^ ^^^^ ^ ^j^) crosses W;q^(v) nOui(w) twice transversely; 

iv) if K > I then the intervals I„ are disjoint. 

Proof. Assertion is immediate from the expression of 77 = $w ° ^'v^w ° ^'v in coordinates: 
ri{x, y) = [x — Klny, y^^ where K = ^Z-,^J^^ > and S = S^S^, > 1 

and assertion In]) follows from Proposition 1^1 

For [m]) let y* be the maximum value of g{x) and let m £ Z such that xq — P^ < ^\ny^ + 2mn. 
We take a — {xo — P^ — 2^1-^)1 K and 6 = (a;o — P^ ~ 2m'K)/K. Then a < b, the second coordinate of 
r]{xo, e'') is less than y*, its first coordinate is P^ + 2rmr, and that of 77(0:0, e") is P^ + 2rmr. Hence the 
curve r?(Xo) goes across the graph of g{x) that corresponds to W^{v) as in figure SI The first coordinates 
of the end points of rjiXn) for the other intervals are PjJ, + 2(?7i + n)TT and P^ + 2(771 + n)Tr and their 
second coordinates are also less than y^,, so each curve r]{Xn) also crosses the graph of g{x) transversely. 

If A' > 1, and since P^ — P^ < 27r < 2Ktv, then xo — P^ — 2rmT — 2Ktt < xa — P^ — 2rmT and hence 
b-2{n + l)n/K < a - 2mT/K, implying that J„ n I„+i = 0. □ 

We are interested in the images of rectangles in /n(v) under iteration by the first return map to 
In{v)\Wil^{v) given by C = *w^v o V- 
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Figure 4. Proof of [mil in Proposition 1101 the transition map rj sends the intervals In into 
curves that cross W^;oc(v) transversely. 



Consider a rectangle i?v = [^'v — '''7 ^'v + '''] ^ 0)^1)] ^ /?i(v) with r > small and < y < y < 1. 
This is mapped by rj into a strip whose boundary consists of two horizontal line segments and two pieces 
of helices. A calculation similar to that in Proposition [lU] shows that there are many choices of y and 
y for which the strip crosses Wfg^{v) transversely near P^. This strip is then mapped by ^w->v into 
In{v) crossing Wfg^{v) transversely. The final strip ({Rv) = ^'w->v ° vi^v) remains close to Wl^^{w). 
Hence the effect of ^ is to stretch i?v in the vertical direction and map it with the stretched direction 
approximately parallel to W^;^^(w). Since Rv has been chosen to contain a piece of W;"^(w), then C(-Rv) 
will cross i?v Repeating this for successive disjoint intervals [y,y] gives rise to horseshoes. 

This is the idea of the proof of the local results of Theorem [3j in a neighbourhood of the connection 
point (i-*y,0), one finds infinitely many disjoint rectangles, each one containing a Cantor set of points 
whose orbits under ( remain in the Cantor set, and hence return to a neighbourhood of {P^,0) for all 
future iterations. In the gaps between these rectangles one finds another set of disjoint rectangles that are 
first mapped by ( into a neighbourhood of the other connection point {P^,0). Repeating the construction 
near the second connection one obtains the switching of Theorem [5] 

It can also be shown that the first return map ^ is uniformly hyperbolic at the points in the rectangle 
Rv, with the choices of y,y above. This means that at each of these points there is a well defined 
contraction direction and this is the main tool in the proof of Theorem 

Since all this takes place inside a positively invariant neighbourhood V'^ , it would be natural to try to 
extend this reasoning to larger rectangles in In{v). We finish this section explaining why this fails. 

Consider a rectangle Rv as above, containing points of the Cantor set. Since 7] expands vertical lines, 
the local unstable manifolds of these points forms a lamination on Rv whose sheets are approximately 
vertical. Now, take a larger rectangle Rv with y S [y, y], y > so as to have the maximum of the curve 
y = g{x) lying in r]{Rv), and increasing t if necessary. The sheets of the lamination still follow the vertical 
direction in the enlarged rectangle, and their image by rj is approximately a helix on Out{w). Changing 
the value of the bifurcation parameter A moves the graph y = g{x) (the stable manifold of v) but does 
not affect the map 77. Hence, by varying A we can get a sheet of the lamination tangent to y = di^), 
say, W^{xo,yo) for some point (xo,yo) in the Cantor set as in figure [S] This breaks the hyperbolicity, 
since it means that W"(a;o,yo) is tangent to W^{v). This phenomenon has been studied by Gonchenko 
et al. [12] — it corresponds to a decrease in topological entropy as in figure [HI As the images of the 
rectangles move down, each time one of them crosses a rectangle a sequence of saddle-node bifurcations 
starts, together with a period-doubling cascade, as on the right hand side of figure [6l 

A more rigorous construction will be made in the next section. 
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Figure 5. Near the point of maximum height of Wi'^^{'v) n Out{w) the unstable manifold of 
some point is tangent to W°(v). By varying A this point may be taken to be in the Cantor set 
of points that remain near the cycle. 




Figure 6. When A decreases, the Cantor set of points of the horseshoes that remain near 
the cycle is losing topological entropy, as the set loses hyperbolicity. This happens when the 
unstable manifold of some point in the Cantor set is tangent to W^(v). 



6. Heteroclinic Tangency 

In this section we show how to find values of the bifurcation parameter A for which W" (w) is tangent 
to Wiv). 

Proof of Theorem^ With the notation of section [SJ the two points (P.^, 0) and (P^, 0) divide the closed 
curve y = g{x) where Wf^^{'v) intersects Out{w) in two components, corresponding to different signs of 
the second coordinate. With the conventions of section [SI we get g{x) > for a; G (^P^,P^y Then the 
region in Out{w) delimited by Wf^^{v) and Wi^^{w) between P^ and P^ gets mapped by Vl/w-i-v into 
the y < component of /n(v), while all other points in Out{w) with positive second coordinates, are 
mapped ito the y > component of /n(v) as in figure [71 The maximum value of the coordinate y for the 
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Figure 7. The transition map ^Pw^v sends the shaded area in Oiit(w), between W;oc(v) 
and W;"(,(w) and with positive second coordinate, into the shaded area in 7n(v) with negative 
second coordinate. Conventions: black hne stands for VK;oj,(w), grey Une is W;o;.(v). The rest 
of the upper part of Out{w) is mapped into the upper part of In{v). 



curve VF;^^(v) n Out{w) is of the order of A, attained at some point {x, y) ~ (a;,„, A) with < x,n < 
(mod 27r). 

Consider now the closed curve where W;"^(w) intersects In{v). For small values of A this is ap- 
proximately an ellipse, crossing Wf^^{v) at the two points {P^,0) and {P^,0), sec figure [31 With the 
conventions of section [3 this is the graph y — h{x) with h{x) > for x e (^P^,P^). In particular, the 
portion of this curve that lies in the upper half of /n(v), parametrised by (x, y), y > 0, may be written as 
the union of two segments cti and (T2 that meet at the point where the coordinate y attains its maximum 
value on the curve. Without loss of generality, let (a;*. A) be the coordinates in In(v) of this point, whith 
P^ <x^< P^ (mod 27r). 

By Proposition [9l the image of each segment aj by i] is a helix on Oui(w) accumulating on W;^^(w). 
Hence, the curve 7/(cri U (T2) is a double helix. The projection of this curve into W^^{w) is regular at all 
points, except for a fold at r]{x^,X) — (x* — K\nX,X^) = {x{X),y{X)), as in figure[8] As A decreases to 
zero, the first coordinate a;(A) of 77(2;*, A) tends to infinity, hence the point 77(2;*, A) makes infinitely many 
turns around the cylinder Out(w). 

On the other hand, y{X) = X^ with (5 > 1, so y{X) decreases to zero faster than A, the maximum 
height of the curve W[^^{v) n Out{w). Therefore, given any small Aq > 0, there exists a positive Ai < Aq 
such that a;(Ai) = Xm and moreover 77(0;*, Ai) — (a;(Ai), ?;(Ai)) lies in the region in Out(w) between 
Wfg^{v) and Wl^^{w) that gets mapped into the lower part of /n(v). For A = Ai, the points on the curve 
W^;"^(w) n /n(v) close to (x*, Ai) are also mapped by 77 into the lower half of Out{w). 

Furthermore, there exists a positive A2 < Ai such that P^ < x{X2) < P^ (mod 2n) and hence 
77(2;*, A2) = (x(A2), ?/(A2)) is mapped by \l/w^v into the upper part of /n(v). Again, for A = A2, points 
on the curve W^^{w) Ci In{v) close to (a;*, A) return to the upper part of Out{w). 

Therefore, for some A*, with A2 < A* < Ai, the image of the curve W;"j,(w) n /7i(v) by the first return 
map to /n(v) is tangent to Wf^^{-v) n In{v), given in local coordinates hy y = 0. This completes the 
proof of Theorem [7] — given any Aq > we have found a positive A* < Aq such that M^"(w) is tangent 
to l^"'(v) for /a,. □ 

7. Conclusion 

For the present study, we have used the symmetry 71 and its flow-invariant fixed-point subspace to 
ensure the persistence of the connections [v w] of one-dimensional manifolds. The symmetry is not 
essential for our exposition but its existence makes it more natural. In particular, the networks we 
describe are persistent within the class of differential equations with the prescribed symmetry. 

As a global structure, the transition from in the dynamics from x = fo{x) to x — f\{x), A w is 
intriguing and has not always attracted appropriate attention. There is a neighbourhood V° with positive 
Lebesgue measure that is positively invariant for the flow of i = f\{x). For A = all trajectories approach 
the network TP. For A > and for a sufficiently small tubular neighbourhood C of any of the 
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Figure 8. In any open neighbourhood of /o in the space of (71 )-equi variant vector fields 
of class on S^, there is a sequence of vector fields />. accumulating of /o whose flow 
has a heteroclinic tangency between Vl^"(w) and W^'*(v). 



Bykov cycles in S'*', almost all trajectories might return to but they do not necessarily remain there 
for all future time. Trajectories that remain in for all future time form an infinite set of suspended 
horseshoes with zero Lebesgue measure. 

For a fixed A > 0, when we take a larger tubular neighbourhood of a Bykov cycle F, the suspended 
horseshoes lose hyperbolicity. While small symmetry-breaking terms generically destroy the attracting 
cycle S", there will still be an attractor lying close to the original cycle. This is the main point of scction[5j 
when local invariant manifolds are extended, they develop tangencies, which explain the attractivity. The 
existence of primary heteroclinic tangencies is proved in section [6] 

Heteroclinic tangencies give rise to attracting periodic trajectories of large periods and small basins 
of attraction, appearing in large (possibly infinite) numbers. Heteroclinic tangencies also create new 
tangencies near them in phase space and for nearby parameter values. We know very little about the 
geometry of these strange attractors, we also do not know the size and the shape of their basins of 
attraction. 

When A — > 0, the infinite number of periodic sinks lying close to the network of Bykov cycles will 
approach the ghost of YP . For the parameter values where we observe heteroclinic tangencies, each 
Bykov cycle possesses infinitely many sinks whose basins of attractions have positive three-dimensional 
Lesbesgue measure. The attractors must lie near W^{'v) U iy"(w) and they collapse into S° as A — )• 0. 
A lot more needs to be done before the subject is well understood. 

Acknowledgements: The authors would like to thank Maria Carvalho for helpful discussions. 
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